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Abstract 

Modified Rayleigh conjecture (MRC) in scattering theory was proposed and 
justified by the author (J.Phys A, 35 (2002), L357-L361). MRC allows one to 
develop efficient numerical algorithms for solving boundary-value problems. It 
gives an error estimate for solutions. In this paper the MRC is formulated and 
proved for static problems. 

1 Introduction 

Consider a bounded domain Del™, n = 3 with a boundary S. The exterior domain is 
D' = M. 3 \D. Assume that S is Lipschitz. Let S 2 denotes the unit sphere in M 3 . Consider 
the problem: 

W 2 v = in D', v = f on S, (1.1) 

v:=0(-) r:=\x\-KX. (1.2) 
r 

Let - := a G S 2 . Denote by Ygipt) the orthonormal spherical harmonics, Yt = Yi m , —£ < 

m < i. Let hi := Y *i+} , I > 0, be harmonic functions in D'. Let the ball Br := {x : \x\ < 
R} contain D. 

In the region r > R the solution to (1.1) - (1-2) is: 

oo 

v(x) = y^cthe, r > R, (1.3) 

the summation in (1.3) and below includes summation with respect to m, — £ < m < £, 
and C£ are some coefficients determined by /. 
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The series (1.3) in general does not converge up to the boundary S. Our aim is to give 
a formulation of an analog of the Modified Rayleigh Conjecture (MRC) from [1], which 
can be used in numerical solution of the boundary- value problems. The author hopes that 
the MRC method for static problems can be used as a basis for an efficient numerical 
algorithm for solving boundary-value problems for Laplace equations in domains with 
complicated boundaries. In [4] such an algorithm was developed on the basis of MRC for 
solving boundary-value problems for the Helmholtz equation. Although the boundary 
integral equation methods and finite elements methods are widely and successfully used 
for solving these problems, the method, based on MRC, proved to be competitive and 
often superior to the currently used methods. 

We discuss the Dirichlet condition but a similar argument is applicable to the Neu- 
mann and Robin boundary conditions. Boundary- value problems and scattering problems 
in rough domains were studied in [3]. 

Let us present the basic results on which the MRC method is based. 

Fix e > 0, an arbitrary small number. 

Lemma 1.1. There exist L = L(e) and q = Q(e) such that 

L(e) 

Il5>(e)/i,-/|| i2(5) <e. ( L4 ) 

£=0 

If (1.4) and the boundary condition (1.1) hold, then 

L(e) 

\\v e ~ v\\ L 2 (s) < e, v e :=^2c e (e)h e . (1.5) 

e=o 

Lemma 1.2. If (1.4) holds then 

\\ Ve - v \\ = O(e) e^O, (1.6) 

where \ \ ■ || := || • \\h™ c (D') + II ' I \l 2 (D';(i+\x\)-i), "f > 1, m > is an arbitrary integer, and 
H m is the Sobolev space. 

In particular, (1.6) implies 

\\v e -v\\ L2{SR) =0(e) e^O. (1.7) 

Let us formulate an analog of the Modified Rayleigh Conjecture (MRC): 
Theorem 1 (MRC): For an arbitrary small e > there exist L(e) and Q(e),0 < 
£ < L(e), such that (1-4) and (1.6) hold. 

Theorem 1 follows from Lemmas 1.1 and 1.2. 

For the Neumann boundary condition one minimizes 1 1 9 ^ e =^ Ce ^ e ^ _ f\\ L2 ^ with re- 
spect to Cf. Analogs of Lemmas 1.1-1.2 are valid and their proofs are essentially the 
same. 

If the boundary data / e C(S), then one can use C(S)— norm in (1.4)-(1.7), and an 
analog of Theorem 1 then follows immediately from the maximum principle. 

In Section 2 we discuss the usage of MRC in solving boundary-value problems. In 
Section 3 proofs are given. 
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2 Solving boundary- value problems by MRC. 



To solve problem (1.1)-(1.2) using MRC, fix a small e > and find L(e) and ce(e) such that 
(1.4) holds. This is possible by Lemma 1.1 and can be done numerically by minimizing 

|| J2o c ^ ~~ /||l 2 (5) := 0(ci, , Cl). If the minimum of <p is larger than e, then increase 

L and repeat the minimization. Lemma 1.1 guarantees the existence of such L and q 
that the minimum is less than e. Choose the smallest L for which this happens and 
define v e := Yle=o c ^- Then, by Lemma 1.2, v e is the approximate solution to problem 
(1.1)-(1.2) with the accuracy 0(e) in the norm || ■ ||. 



3 Proofs. 

Proof of Lemma 1.1. We start with the claim: 

Claim: the restrictions of harmonic functions hi on S form a total set in L 2 (S). 

Lemma 1.1 follows from this claim. Let us prove the claim. Assume the contrary. 
Then there is a function g ^ such that f s g(s)hi(s)ds = W > 0. This implies 
V{x) := f s g(s)\x — s| — 1 cZs = Wx G D'. Thus V = on S, and since AV = in D, one 
concludes that V = in D. Thus g = by the jump formula for the normal derivatives of 
the simple layer potential V. This contradiction proves the claim. Lemma 1.1 is proved. 
□ 

Proof of Lemma 1.2. By Green's formula one has 

w e (x) = / w e (s)Gn(x, s)ds, \\ui e \\ L 2(s) < e i w e :=v e — v. (3.1) 
Js 

Here N is the unit normal to S, pointing into D', and G is the Dirichlet Green's function 
of the Laplacian in D': 

V 2 G = -S(x - y) in D', G = on S, (3.2) 

G = 0(-), r^oo. (3.3) 
r 

From (3.1) one gets (1.7) and (1.6) with H^ c (D')—noYm immediately by the Cauchy 
inequality. Estimate (1.6) in the region B' R := M 3 \ Br follows from the estimate 

\G N (x,s)\<^— \x\>R. (3.4) 

In the region Br\D estimate (1.6) follows from local elliptic estimates for w t := v e — v, 
which imply that 

\\we\\L2(B R \D) < ce. (3.5) 

Let us recall the elliptic estimate we have used. Let D' R := Br\D and Sr be the boundary 
of Br. Recall the elliptic estimate for the solution to homogeneous Laplace equation in 
D' R ( see 0. p.189): 

< c[\\w € \\ L 2 {Sr) + \ \w e \\ L 2 [s) ]. (3.6) 
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The estimates ||tu e | |l2( 5h ) = 0(e), | |^y e | |x, 2 (S) — 0(e), and (3.6) yield (1.6). Lemma 1.2 is 
proved. □ 
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